Non-adiabatic geometric quantum computation with trapped ions 
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We propose a non-adiabatic scheme for geometric quantum computation with trapped ions. By 
making use of the Aharonov-Anandan phase, the proposed scheme not only preserves the globally 
geometric nature in quantum computation, but also provides the advantage of non-adiabaticity that 
overcomes the problem of slow evolution in the existing adiabatic schemes. Moreover, the present 
scheme requires only two atomic levels in each ion, making it an appealing candidate for quantum 
computation. 



Conventionally, the controllable operations in quantum 
computation (QC) are achieved on the basis of dynamical 
origins of quantum state evolutions pj. In recent years, 
for the purpose of being fault-tolerant to certain types of 
computational errors, there are considerable interest in 
exploiting the possibility of implementing quantum com- 
putation by geometrical means, which have been termed 
as holonomic quantum computation 0, 0, 3 • Depending 
on the degenerate property of the eigenspace of the gov- 
erning Hamiltonian, the holonomy can be cither a simple 
Abelian Berry phase factor Q or a general non-Abelian 
unitary transformation It has been shown that the 
universal quantum computation can be implemented in 
principle by holonomies OJEll Up to date, several ex- 
perimental proposals have been suggested for geometric 
quantum computation, using such as the nuclear mag- 
netic resonance 0, superconducting nanocircuits 0, |9| , 
trapped ions [Tol lll|. and nonlinear optics |T^. 

The principle of the aforementioned geometric QC is 
rooted in the adiabatic evolution of quantum system, 
which may thus imply a slow computing speed. The adi- 
abatically slow evolution may also challenge the sustain- 
ment of required coherence in QC. Therefore, geometric 
QC based on non-adiabatic evolution should be desirable. 
Very recently, Wang and Keiji suggested to exploit a non- 
adiabatic evolution to realize geometric QC in NMR sys- 
tem and superconductor nanocircuits [l3|. Indeed, ge- 
ometric phase exists in non-adiabatic evolving quantum 
systems, which is in fact the Aharonov-Anandan (A- A) 
phase [14J. Strictly speaking, the A- A phase depends 
on certain dynamical quantities such as the rotating an- 
gular s pee d of external (magnetic) field or state vector 
[l&lril lTsI]. In this sense, the A- A phase differs from the 
adiabatic Berry phase. However, the dependence of the 
A-A phase on the angular speed is through the closed 
path loop depicted by the ending point of state vector, 
and hence is global in nature that largely retains the ge- 
ometric sense of A-A phase. Accordingly, quantum com- 
putation based on the non-adiabatic A-A phase can be 
reasonably regarded as a kind of geometric QC. 

In this work we propose a scheme for non-adiabatic 
geometric quantum computation with trapped ions. Be- 



sides removing the drawback of the slow adiabatic evo- 
lution, the proposed non-adiabatic scheme holds addi- 
tional merits as follows. Firstly, there is no need to 
design the reverse evolving path to eliminate dynamical 
phases that occur in the Berry phase-based, adiabatic ge- 
ometric QC operations in nondegenerated systems jj], |8| • 
Secondly, in comparison to the existing fully holonomic 
QC schemes !). fjjf . the present one does not involve the 
complicated construction of the degenerate eigenspace of 
driving Hamiltonian. The above two merits stem from 
the following observations. In a non-adiabatic quantum 
evolution, the geometric A-A phase is in general accom- 
panied by a dynamical phase. However, if the evolving 
path is such designed that along it the state vector is 
always perpendicular to the driving (magnetic) field, the 
resulting phase factor after a non-trivial cyclic evolution 
will be purely geometric. This feature has been exploited 
by Suter et al in their seminal experiment for demonstrat- 
ing the A-A phase Finally, only two atomic levels 
of each ion are needed in our scheme. This merit alone 
is attractive, since the originally proposed ion-trap QC 
scheme required three levels |19| and the recently pro- 
posed holonomic ion-trap QC required four levels plfll . flll | . 
We notice that valuable efforts on improving ion-trap QC 
protocol by using only two levels have been carried out 
in dynamic schemes Eel l2lt 12^ ]. In particular, the tech- 
nique proposed in Ref. I22L which effectively couples the 
electronic states of a pair of ions by virtually exchanging 
phonons, is similar to our present one. The major contri- 
bution of our work is to perform QC by geometric means. 
Viewing that the work of Ref. |22f has in fact extended the 
ion-trap QC scheme to finite temperature, our geomet- 
ric scheme may also hold to similar regime, although in 
the following we would restrict our description at zero 
temperature limit. 

Model Description. For quantum logic with trapped 
ions, we assume that each ion has two relevant internal 
states |0) and 1 1) with energy separation loq, and, as usual 
[l9|. can be selectively addressed by lasers. Consider, for 
instance, the jth ion being exposed to a traveling-wave 
laser field E(z) = E cos(k • z — ui^t + <P) with frequency 
u>l, wave vector k, and phase cf>. Here z = zoi(a + a^) is 
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FIG. 1: Schematic diagram for geometric rotation of state 
vector around the effective magnetic field. Since the state is 
always perpendicular to the field, there is no dynamical phase 
accumulation during the evolution. 



A-A phase under the following operations: (i) Switch- 
ing on a 7r-pulse with laser phase = 0, the state 
|+) rotates around an effective magnetic field fii = 
{ujj,0,(ujQ — ul)/2}, from |+) in the e y direction to 
|— ) in the — e y direction along the curve ACB on the 
Bloch sphere; see Fig. 1. (ii) Suddenly changing the laser 
phase to = 7r, after another 7r-pulse, the state |— ) ro- 
tates back to |+) around fi 2 = { — w^O, (ujq — wl)/2} 
along curve BDA on the Bloch sphere. According to 
the A-A phase theory, after the above cyclic evolution, 
the state |+) will acquire a geometric phase e 27 , with 
7 = 4 arctan[2Wj/ (cjq —&>£)]• Note that during the above 
operation, the state vector keeps always perpendicular 
to the effective magnetic field, thus no dynamical phase 
is accumulated in the evolution. Similarly, the state |— ) 
will acquire A-A phase e~ n at the same time. 

Now consider the evolution of logic states |0) = 
^§(|+) - |-)), and |1) = -^(|+) + |-)). After the above 
operations, they evolve to states 



the center of mass coordinate of the ion in terms of the 
phonon raising (lowering) operator a^(a) and zero-point 
spread zq = (h/2Muj) 1 ^ 2 , where M is the total mass of 
the ion chain, and u> the phonon frequency. The resulting 
Hamiltonian reads 



H.c. 



(1) 



Here, the atomic Pauli operators <j z = |1)(1| - 
(j + = |1)(0|, and cr~ = |0)(1|, are introduced. In 
Eq. 0} ujj is the Rabi frequency, and rj = (k • z)zq is the 
Lamb-Dicke parameter, which accounts for the coupling 
strength between internal and motional states. In our 
geometric scheme, the phonon plays no roles in one-bit 
operation, but it does involve in two-bit logic gate imple- 
mentation as it effectively couples two ions together. 

One-Bit Rotation via Non-adiabatic A-A Phase. As 
just mentioned, in single bit logic operation, we do not 
need the phonon- assisted dynamical flipping. This can be 
the case under situations such as (i) at zero temperature 
(no phonon excitation), if ujl < ujq; and (ii) at finite 
temperature, but luq — ^ fiu) (n is an integer). In 
the absence of phonon participation, Eq. can be re- 
expressed as 
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in the rotating frame (with angular velocity w/,e z ). Here 
notations have been introduced for the effective magnetic 



|0) -> cos 7 |0) +sin 7 |l), 
|1) — > cos7|l) — sin7|0). 



(3) 



We see here the geometric A-A phase plays a role of rotat- 
ing the logic states. Particularly, complete state flipping 
between |0) and |1) can take place at 7 = tt/2. Note also 
that the possible value of 7 ranges from to 27r, implying 
the ability of arbitrary rotation between |0) and |1). 

The state evolution of the performed qubit described 
by Eq. J2Jl is expressed in the rotating frame with fre- 
quency lol, in which other free (not performed) qubits 
would have relative phase accumulations in the non- 
resonant case of u>l 7^ luq. Conventionally, a more conve- 
nient choice is to express states in the interaction picture 
with respect to Hq = ^ J2j a j (equivalently, a rotating 
frame with frequency ujq around z-axis). Accordingly, 
Eq. J2J can be recast in the interaction picture as 
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(4) 
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field flj = {ujjcos* 
Pauli operator a , _ , , , ; - . , 

To show how rotation between logic states |0) and |1) 
can be performed by geometric means, we first consider 
the geometric evolution of the eigenstates of cr v , defined 
by cr y \±) = ±|±), under appropriately designed laser 
sequences. Specifically, the state |+), for instance, will 
complete a cyclic evolution and acquire a non-adiabatic 



where ljd = ^0 ~ and t is the total operation time on 
the performed qubit. 

Interestingly, the above state rotation (flipping) is per- 
formed by non-resonant pulses via geometrical means. 
Now we show that by resonant pulses, a pure phase 
shift gate of single qubit can be geometrically real- 
ized. In resonant case, the laser-frequency-associated 
rotating frame coincides with the interaction picture 
defined by 7io, in which the (jth) qubit Hamiltonian 



sin^, (uiq — cj.l)/2}, and the vector reads H^' = Cuj [a^e 1 ^ + H.c] = Slj ■ &j, with flj = 



{ujjCOS(f>, ujj sum/), 0}. We see that in the rotating frame the 
effective magnetic field constantly has zero z-axis com- 
ponent for arbitrary laser phase 0, i.e., it lies in the x—y 
plane. To realize the single bit phase gate, we first switch 
on a 7r-pulse with laser phase at certain value, say, —cf>o- 
The logic state |0) and |1) would rotate around the effec- 
tive magnetic field {cl^cos^oj — Wjsin^o, 0} to |1) and |0), 
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respectively. Then, suddenly changing the laser phase to 
<fio, after another 7r-pulse, the logic states rotate around 
the effective magnetic field {o^cos^o, QjSm<f>o, 0}, and re- 
turn to the original states |0) and Associated with 
this two step cyclic evolution, the logic states will respec- 
tively acquire geometric A-A phases as 
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where 7 = 40o- With the help of this phase shift oper- 
ation, the additional phase factor in Eq. J3J can be can- 
celed out by properly choosing <pQ. More importantly, 
together with this phase shift gate, the qubit state rota- 
tion Eq. (0J constitutes a complete logic set for arbitrary 
single qubit operation. 

Two-Bit Gate via Non-adiabatic A-A Phase. To com- 
plete the universal gate for quantum computation, non- 
trivial two-bit gate such as the controlled NOT, or equiv- 
alently, the conditional phase shift (CPS) gate, would be 
required. Below we show how the CPS gate can be imple- 
mented via geometric means. Consider two qubits (e.g. 
the jth and kth ones) irradiated by two lasers with fre- 
quencies and u>l,2, and phases <f>\ and <p2- By setting 
(jJl.i > 1^0 an d wl,2 < 1^0 > an d correspondingly denoting 
the detunings by <5i = (jJls — ^0 and 82 — ~ &L,2, the 
effective coupling between the two-bit states |00) and |11) 
can be established via virtually exchanging phonons, and 
the resulting two-bit effective Hamiltonian reads [Z| 
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Here notations {|1) = |11), |2) = |00),|3) = |10),|4) = 
|01)} are introduced for the two-bit computational basis 
states. Up to the first-order expansion of the Lamb-Dicke 
parameter 77 in Eq. the effective coupling strength 
can be obtained via second-order perturbation theory as 



9jk = 9j9k{^ 



where 



phonon involved Rabi frequency of single ion transition. 
The four basis-state energies E m (m = 1, • • • , 4) contain 
also the ac Stark shifts. The effective interaction couples 
only between |1) and |2), but leaves |3) and |4) inactive 
with respect to the laser operation in study. This effec- 
tive two state dynamics (e.g. Rabi oscillations) has been 
demonstrated in Ref.l22lbv numerical simulation starting 
from the original Hamiltonian. As a result, for two-bit 
operation, one can focus on state evolution in the sub- 
space {|1), 1 2)} in which the two-bit Hamiltonian (6) is 
recast to a similar form as the single qubit |23j 



j(k)V is the one- 
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where $ = 4>i + 4>2, an d the two-bit Pauli matrices 
are introduced _as E z |l(2)) = ±|1(2)), S+|2) = |1), 
and E~|l) = 1 2). Note also that the Hamiltonian Q 



has been expressed in the rotating frame with respect 
to Y7' with the rotation frequency ujla + ^l,2) thus 
ujd = {E\ — E2) — + u>l,2)- Simple comparison 

of Eq. |(7J| with Eq. J2J) indicates that an arbitrary rota- 
tion between |1) and |2) can be performed geometrically 
as that in the single bit case. 

In the spirit of pulse-sequence operations in dynamic 
scheme based on an XY spin model |2J, |25|, |26j , we shall 
in the following show that the CPS gate can be imple- 
mented geometrically as 



CPS 



^i-K/A^i-KriyCr j/3 iirn.k-(T k /3 —iircr% /2 
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Here the vector Pauli operator er = (a x ,a v ,a z ), unit 
vector iij = (1,1,-1)/V3, and n fe = (1, -1, l)/\/3. In 
the interaction picture with respect to the original free 
ion Hamiltonian, the two-bit rotation operator Ujk{cx) is 
defined via U jk (a)\01) = |01), U jk (a)\10) = |10), and 
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U jk {a)\QQ) = cos (|) |00) + isin 

U jk (a)\ll) = cos (f) |11) + zsin (|) |00). (9) 

Since the arbitrary rotation and phase shift gate of sin- 
gle qubit have been built up, the single bit operations 
in Eq. © can be readily implemented via geometric 
means by properly combining the single bit logic ele- 
ments. Viewing the similarity between the two-bit re- 
duced Hamiltonian (JZJ) in the subspace {|00), |11)} and 
the single bit Hamiltonian @, the two-bit rotation of 
Eq. can be straightforwardly implemented by the fol- 
lowing two-step procedures: 

(i) In the two-bit rotating frame with frequency u>£ 1 + 
u>l,2 around E z , performing cyclic evolution for the eigen- 
states of E y by controlling the laser phases $ similarly 
as in the single bit case, one can geometrically rotate the 
states 1 00) and |11). Expressed in the interaction picture 
which also corresponds to H = £?i|ll)(ll| + £ 2 |00)(00|, 
this operation can realize the following state transforma- 
tion 

|00) -> e _ii>DT/2 cosr|00) + e^ T / 2 sinr|ll), 

cosr|ll)-e- iliDT / 2 smr|00). (10) 
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Here Y is the geometric A-A phase determined by the 
evolution contour of the two bit state vector, and (ijjT is 
the detuning-induced phase accumulation. 

(ii) Tuning the laser frequencies in resonance with the 
two ions |23], i.e., ujd = 0, a phase-shift gate associat- 
ing with 1 00) and |11) can be implemented to cancel out 
the phase factors in Eq. I|10l) . In the resonance case the 
+^l,2 rotating frame coincides with the interaction 
picture of Ho- As the one- bit resonant case, the effective 
magnetic field corresponding to Eq. (JJJ now lies in the x- 
y plane since Cjd — 0. By successively choosing two dif- 
ferent values of the laser phase $, one can perform two 7r 
rotations on the states |11) and 1 00) around the effective 
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magnetic fields, and readily generate the A-A geometric 
phases, e~ lT and e %T , for states 1 11} and 1 00) , respectively. 
Now, after a phase-shift operation with T = ludt/2 — 7r/4, 
Eq. l|lUfl becomes 

|00) -► e~ l7r/4 [cosr|00) +isinr|ll)] 

|11) -> e™ /4 |cosr|ll) +isinr|00)| . (11) 

This is identical to Eq. (J5J) , except for the additional 
global phases. Obviously, these global phases will not 
appear if we first generate a phase shift of e m l A on |00) 
and e~ l1T / 4 on |11) at the same time, by the phase-shift 
gate just described, prior to the operation of Eq. H10JI . 

We have thus realized the two-bit gate Ujk as defined 
in Eq. @ . Together with the arbitrary one-bit operations 
(rotation and phase shift) described earlier, we can now 
readily implement the important CPS gate Eq. JHJ, whose 
role is to transform 1 1 1) — * e I7r |ll), while to keep other 
computational two-bit basis states unchanged. 

Conclusion and Discussion. The proposed non- 
adiabatic geometric QC scheme based on the A-A phases 
is expected to overcome several drawbacks of the adia- 
batic schemes, namely, the slow evolution, need of refo- 
cusing to eliminate the dynamical phases, and continuous 
control over many fields to construct non-trivial loops in 
the parameter space. Viewing that the trapped-ion is 
one of the best exploited systems for quantum computa- 
tion, and that the proposed scheme requires a relatively 
simple atomic level configuration, we suggest, as a first 
step, to exploit it as an interferometer for principle-proof 
of the non-adiabatic geometric A-A phase discussed in 
this work. We believe that the interference associated 



with the non-adiabatic A-A phase can be readily demon- 
strated by experiment in ion-trap systems. 

As a possible QC architecture, the elementary op- 
eration steps in the proposed non-adiabatic geometric 
scheme are comparable to its dynamic counterparts. 
Specifically, the time scales for both the one-bit and two- 
bit geometric operations are about the same as those in 
the dynamic operations. Concerning the possible fault- 
tolerance, in the adiabatic case, quantum logic is toler- 
ant to certain types of errors such as the field fluctua- 
tions that preserve the loop area in parameter space, i.e., 
the Berry phase. Similarly, in the non-adiabatic case, 
the A-A phase is of error tolerance to any fluctuations 
around the state-evolving path that preserve the path 
loop area. In principle, there exist many possible driv- 
ing field deviations that can preserve the state path loop 
area. However, the most natural and possible errors ap- 
pear to be random (but small) fluctuations of the laser 
phase, frequency, and coupling strength to the atomic 
levels, which equivalently result in fluctuations of the ef- 
fective magnetic field. The global A-A phase is expected 
to be largely immune from this type of errors, and, at 
the same time, the net dynamic phase accumulation is 
approximately zero due to the cancellation of the posi- 
tive and negative contributions. 
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